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1. Introduction. 



The purpose of this paper is to present simple models of the flow of 
officers through P-coded billets in order to determine yearly flows into 
graduate education to meet future billet requirements. Recent significant 
reductions in P-code billets in certain disciplines have led to serious prob- 
lems in managing current inventories, in determining future educational inputs, 
and in determining how educational institutions should adjOst to severely 
reduced inputs in a short period of time. 

This report is written in six sections of which this is the first. In 
sections 2 through 5 it is assumed that the requirements for officers with 
graduate degrees is known (in section 6 some comments are made regarding 
determination of these requirements). Section 2 describes a steady state flow 
model in detail. An interactive computer program based on this model is 
described in Appendix (i). Section 3 describes two models which deal with 
the transient problems created with billet reductions. The first model leads 
to unsatisfactory cyclic quotas which result from an unrealistic (but not 
obvious) assumption in the model. In the second model we modify this hidden 
assumption and obtain much smoother quotas. The mathematical formulation of 
these models are given in the Appendix in sections (ii) and (iii) together 
with details of APL functions based on the models. 

Section 4 contains comments and reasons for smoothing beyond the narrow 
mathematical considerations of section 3. Some of the factors discussed are 
difficult to quantify and so little mathematical reasoning is given. However, 
further smoothing of the quota from past and current levels to the new steady 
state levels given by the model in section 1 may be very important to the 
quality of graduate education. 



Section 5 discusses problems in using the current inventory/billet 
ratio for a given P-code as an indicator of the "health" of that P-code 
community. Modifications of the current method are suggested and examples 
are given to show how the current ratio can be misleading in a transient 
stage caused by billet reduction. 

Finally, in section 6, some comments are made concerning P-code billet 
requirements. These requirements drive the whole quota system and the models 
of sections 1 through 3 are for nought unless the billet requirements given as 
input to the models truly represent the requirement. Although it is not the 
purpose of this report to determine requirements it would not be complete 
without mentioning some of the problems associated with this difficult subject. 
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2 . The Steady State Model . 



The model developed in this section equates input rates and output 
rates to determine the basic flows through the officer billets. Minimum 
steady state flow rates (officers per year) into graduate education are deter- 
mined which will meet future billet requirements. 

In order to simplify notation, consider the P-code and designator fixed. 
The subscript i or j on the variables indicates rank, as follows: 



Subscript Rank 

1 LT 

2 LCDR 

3 CDR 

4 CAPT 



Let be the number of P-coded billets in rank i (for the fixed 

P-code and designator), and T^ the tour length in these billets. If the 

B^. billets have varying tour lengths, then T^. can be thought of as an 
average. 

Current career paths call for a time lag from the time an officer 
graduates to the time he serves in a P-coded billet. Many different paths 
are possible in this period and much turbulance can take place. This turbu- 
lance is summarized in a single parameter for each P-code and designator: 

= Fraction of officers who enter graduate education to meet a future 

4 billet requirement in rank i, who are still in the Navy and 

/ eligible to meet that requirement when it occurs, i = 1,2, 3, 4. 

('(ib vf r \ f j jo tjf - 

Note that the officer student may not be in rank i when he enters 
school. Typically his rank would be (i-1) with the current intervening 2-3 
year tour between graduation and entrance into a P-coded billet. Clearly 
is a function of policy as well as rank, P-code and designator. If a student 
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immediately entered a P-code billet on graduation then a^. would be 1. If 
a student entered school as a LT to meet a P-code billet requirement as LCDR 
y ""in five years (th^current BuPers Model assumes this), and selection from LT 
to LCDR was 0.75, then would be closer to 0.75. 

Because of the complexities and interactions between the promotion 
structure, career paths, and rotation dates, it is possible that not all P- 
coded officers will get to serve a P-coded tour in rank. Again, the numerous 
alternative paths which officers can take cannot (and should not) all be 
treated as separate possibilities; rather they are summarized by a single 
parameter; ,f 

3- = Fraction of those available to serve a P-code tour in rank i 
who get to serve such a tour. 

In a perfect system 3.j would be 1. In practice it would be somewhat 
less than 1. Later we demonstrate how the parameters and 3.j can be 

used to give bounds on the student input flow. 

For ranks above LT it is possible (and desirable) that officers serve 
second or third tours in P-coded billets. We let 



Y.j = Fraction of those eligible to serve in a P-code billet in rank i 
who are available to serve a P-code tour in rank i + 1. 



The network of flows and inventories is shown in Figure 1. The basic 
, flow rates for each rank are determined by the billets and tour lengths. Thus 

J . 

] if in rank i there are B. billets with an average tour length T., then 

^i ^ 

of these billets become vacant each year (assuming steady state flows). 

Thus we must have an input flow equal to this number to meet the billet require- 
ments. 
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Figure 1 : Network of Steady State Flows Through Officer P-Code System . 
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Let X^. be the input per year into graduate education necessary to 
meet the billet requirements in rank i, i = 1,2, 3, 4. Starting with i = 1 
it is easy to see that the flow into graduate school to meet the LT billet 
requirements, is given by 



X, = 



1 



ft , 



(1) 



The reader can check that with this input, if remain in the system 

to fill a billet, and a fraction of these get to serve in the billet, 

^1 

then the flow into the LT billets is the outflow from these billets. 

‘l 

We move now to i = 2, the LCDR billets. Input into these must be 
B B 

Y~ per year. Now w v officers are in service at the end of the LT P-code 
*2 



B,T, 



tour. Of these a fraction y-| a>^e available to serve in LCDR billets some 



time later. Now if 



y,b, 



B 



o ' Y ^ then not all these get to serve a second tour 
'^I'l *2 



and no new input to graduate school is required for future LCDR billets. How- 



ever, if 



Y,B 



B. 



ri . "2 



e,T, 



j , then we need some entry into graduate school to make up 
*2 



the deficit in the flow rate. It follows that 

1 



'2 = a.,3 



Max 



2^2 



0 , 



®2 



e,TiJ 



or, by substituting (1) in the right-hand-side, 

\ Max(0, - a^y^X^). 

Similar arguments show that 



( 2 ) 



1 



X, = Q 
3 a»3 



Max 



and 



3^3 

1 



^3 

T3 ■ ^“l^l^2^1 “2^2^2^ 



(■64 

r’ T7 ■ ^“l^l^2^3^1 “2^2^3^^2 ^ ^3^3^3\ 



^4 a,3 



Max 



4^4 



(3) 



(4) 



6 



Finally, the total input per year to graduate school for the particular 



P-code and designator is ^2 ^ ^3 ^ ^4’ 

Example : The current (September 1974) billet requirements for code 8510P, 

unrestricted line are given below in Table 1 together with three year tour 
lengths which are assumed for this example. The values of a, 3 and y 
aro chosen assuming three year tours between graduation and entering a P-code 
billet, 75% selection of LT's to LCDR's, 70% selection of LCDR's to CDR's, 
and 50% selection of CDR's to CAPT's. We have also assumed a 90% utilization 
of graduates in their first tour, and 100% re-utilization (this last assumption 
leads to Y.j_i The values of X.j to calculated using (1) to (4) 

are shown in Table 1. 



Rank 

i 


Bi 1 lets 

B. 

1 


Tour Length 

T. 

1 


a . 
1 




^i 




1 


35 


3 


.95 


.90 


.75 


14 


2 


67 


3 


.75 


.90 


.70 


19 


3 


59 


3 


.70 


.90 


.50 


5 


4 


12 


3 


.50 


.90 


— 


0 


Total 


173 


- 


-- 


— 


— 


38 


Table 1 


: Example 


of Quota Determination 


Using 8510P 


URL Data 





The parameters a, 3, Y are of course subject to interpretation and 
are affected by policy changes such as changes in career paths. Suppose we 
wish to find a lower bound on the input to meet the billets in this example. 
This is done by setting a^. = 3.j = Y.j = 1 for all i. What this says is that 
everyone sent to school as a LTjg or above will serve in a P-coded tour at 
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every opportunity and there will be 100% selection all the way to CAPT. With 
these extremes we get = 12, X 2 = 10, X^ = X^ = 0 for a total input of 
22 per year. This lower bound would never be sufficient to meet requirements 
in a real system, but might serve as useful information when trying to estimate 
the effects of uncertainty about the a, 3 and y parameters. Suppose now 
that we keep the a^. and 3^- fixed as in Table 1, but reduce the tour lengths 
to 2.5 years and have no one complete two or more P-coded tours. Thus = 0 
and the new input would be: X^ = 16, X^ = 40, X^ = 37, X^ = 11, for a 

total of 104. 

Clearly both these examples are extremes, but are given to illustrate 
how the model can be used to determine the effects of different policies. 

Using the same data as in Table 1, but with 3.j =1.0 i = 1,2, 3, 4, we 
obtain Table 2. Thus by using all P-coded personnel at every opportunity the 
total quota is reduced from 38 to 36. Suppose in addition to all ^-j “ 1 we 
plan on two 3-year tours in the rank of LCDR. We then set ^2 ~ ^ repeat 
the calculations. The results are shown in Table 3. Note that the total 
quota is decreased to 32, but the "mix" of the four inputs changes considerably. 

An interactive computer program has been written which enables the user 
to change the various parameters for a given P-code and designator at a terminal. 
The program is written in APL and is given in the Appendix together with an 
explanation of how it is used and an example. 
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Rank 

i 


Bi 1 lets 


Tour Length 


“i 




^i 




I 


35 


3 


.95 


1.0 


.75 


12 


2 


67 


3 


.75 


1.0 


.70 


18 


3 


59 


3 


.70 


1.0 


.50 


6 


4 


12 


3 


.50 


1.0 


— 


-- 


Total 


173 


— 


— 


— 


-- 


36 


Table 


2: Quota calculation using 


Table 1 


data, but with 


3,- = 


1. 



Rank 

i 


Bi llets 


Tour Length 

T. 

1 


a. 

1 




^i 




1 


35 


3 


.95 


1.0 


.75 


12 


2 


67 


6 


.75 


1.0 


.70 


3 


3 


59 


3 


.70 


1.0 


.50 


17 


4 


12 


3 


.50 


1.0 


— 


-- 


Total 


173 


-- 


— 




-- 


32 



Table 3 : Quota calculation using Table 2 data, but with = 6 . 



3. Transient Problems and Quota Smoothing . 



The model in section 1 assumes that the number of billets in a given 
rank, P-code, and designator are constant over some reasonable time period. 

The model does not have enough detailed structure to enable us to examine the 
transient effects of a change in the number of billets. Rather large decreases 
in P-coded billets have recently occurred and the effects of such changes on 
the quota are the subject of this section. 

Before giving a mathematical flow formulation of the transient problem 
a simple example is used to illustrate the concepts. Consider a fictitious 
system which, up to and including the current planning period, had 90 billets 
requiring school education. For simplicity we assume that people are sent to 
school in year 1 of their career for a 1 year school program, and immediately 
follow this with three years in one of the 90 billets. We assume that everyone 
graduates and stays in the system at least four years. 

Example 1 : A Simple Illustrative Example . 

This simple example is illustrated in Figure 2. Planning periods 
increase down the page and period 0 is the current period. Since we assume 
we have been running in steady state the current flowrate is 30 people per 
year out of the billets. Thus we have 30 people in school in period 0 ready 
to fill the vacated bilets in period 1. The four numbers in the top line boxes 
are a legacy of previous policy and cannot be changed. 

Assume that a 40% reduction in billets occurs in period 1 and that the 
number of new billets will be constant at 54. Since there are 90 educated 
people in period 1 we have an unavoidable excess inventory of 36 people. Our 
problem now is to calculate the new school inputs in periods 1,2,3,... . 
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Planning 

Period 



0 



1 



2 



3 



4 



5 



6 



7 



School Tour 



1 2 3 4 



30 




30 


30 


30 




0 




30 


30 


30 










24 




0 


30 


30 












30 


24 


0 


30 




0 




30 


24 


0 








24 




0 


30 


24 




30 




24 


0 


30 














30 


24 


0 



Total 
Bi 1 lets 



90 



54 



54 



54 



54 



54 



54 



54 



Excess 

Inventory 



0 



36 



6 



0 



0 



0 



0 



0 



Figure 2 . 



Personnel Flows in a Simple Transient Example . 
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In period 2 the total inventory will be 60 plus the input from those 
trained in school in period 1. Since 60 already exceeds the 54 billets the 
school input in period 1 is zero. In period 3 the total inventory will be 
30 plus the input from school in period 2. Since we need 54 the school input 
in period 2 must be 24. 

It is easy to calculate the inputs for this example. For periods 1 
through 6 they are 0, 24, 30, 0, 24, 30. The reader can see immediately 

that a cyclic input results. This is an extremely undesirable feature for 
planning purposes. Firstly, it adversely affects the school which must try 
to meet widely varying inputs with a relatively stable faculty. Secondly, it 
leads to gross inequities in educational openings between year groups. 

This cyclic input feature is not restricted to our simple model, but 
is a result of the underlying arguments. These arguments hide some unwitting 
and unrealistic assumptions to which we return later. First we demonstrate 
the same cyclic feature of the quota derived for a realistic example. 

Example 2 : A Realistic Example Showing Cyclic Quotas . 

Table 4 gives the basic data used in the example. It is based on real 
data for the 8510 P-code of URL officers and billets in November 1974. The 
current inventory (period 0) is given in column 2 by year group. Columns 3, 

4 and 5 give the expected additions to the inventory in planning years 1, 2 and 
3 from students currently enrolled in the 360 curriculum at Monterey. Column 
6 gives the assumed continuation rates from the year given by the row, to the 
next year. Column 7 gives the total P-coded billets by rank, which are the 
same as those used in the example in section 1. The horizontal lines indicate 
that LT billets are held by officers with 8-10 years commissioned service, 

LCDR billets by officers with 14-16 years, CDR billets by officers with 19-21 
years, and CAPT billets by officers with 24-26 years of service. 
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Years 

Commissioned 
Servi ce 


Current 

Inventory 


Student 


Additions 


Continuation 

Rates 


P-Coded 

Billets 




Yr 1 


Yr 2 




1 










1.00 






2 










1.00 






3 


4 








1.00 






4 


6 


1 


2 


2 


0.98 






5 


6 


1 


4 


4 


1.00 






6 


9 


4 


8 


2 


1.00 






7 


13 


3 


8 




1.00 






8 


21 




8 




1.00 


\ 




9 


36 




3 


2 


0.75 


[ 35 


LT Billets 


10 


35 


4 


3 


3 


1.00 


) 




11 


51 


4 


2 




1.00 






12 


43 




1 




1.00 






13 


42 


1 






1.00 






14 


33 


1 


1 


1 


1.00 


\ 




15 


37 




1 




1.00 


[ 


LCDR Billets 


16 


30 








0.70 


) 




17 


25 








1.00 






18 


25 








1.00 






19 


21 








1.00 


) 




20 


16 








1.00 


1 


CDR Billets 


21 


11 








0.50 






22 


5 








1.00 






23 


15 








1.00 






24 


5 








1.00 






25 


6 








1.00 


( 12 


CAPT Billets 


26 


6 








1.00 


) 




27 


2 








1.00 






28 


2 








1.00 






29 


2 








1.00 






30 


6 








0 






Total 


512 















Table 4: Basic Data Input for Transient Quota Model Calculation. 
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A mathematical formulation and description of the model is given in 
Appendix (ii) together with a computer program listing in APL and a sample 
example of the input and output. The results of the quota calculations for 
15 years using the data in Table 4 are shown in Table 5. We have assumed a 
5-year lag between entering school and entering a P-coded billet. Thus group 
1 input would enter school in their third year of service, group 2 in their 
ninth year, group 3 in their fourteenth year, and group 4 (if any) in their 
nineteenth year. 

A cyclic trend is clearly evident in groups 2 and 3. The total quota 
is also cyclic, and fluctuates from about 20 to 45 after an initial period. 

The reader can check that the averages of the 3-period cylces give 13, 17, 

6, 0 respectively, for the four groups, which agree closely with the results 
of the steady state model (section 1) in Table 2 with 3^. = 1.0, i = 1,2, 3, 4. 

It is implicitly assumed in the transient model that all educated offices 
will serve in a P-coded tour at every opportunity . 

Example 3 : The Simple Example with Different Assumptions . 

We now return to our simple example where the 90 billets are reduced 
to 54. Although Figure 2 is simple and seems to demonstrate the correct flows, 
they are correct only if a hidden assumption is valid. This assumption is 
illustrated in Figure 3. Each box in each planning period now contains two 
numbers. The lower number in each box is the number of educated people filling 
P-coded billets. The upper number is the number of people educated, but not 
holding a P-coded billet. By looking at these two different communities an 
unrealistic assumption which is present, but hidden, in Figure 2 is demonstrated. 

In the current period 0, since there are 90 billets there is no excess 
inventory of people in their second, third or fourth year. In period 1, 36 
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Planning 
Peri od 


School Quotas 


GRP 1 


GRP 2 


GRP 3 


GRP 4 


TOTAL 


1 


7 








7 


2 


12 








12 


3 


12 


22 






32 


4 


13 


7 






20 


5 


14 


5 






19 


6 


12 


33 






45 


7 


13 


12 


17 




43 


8 


13 


6 






20 


9 


12 


33 






45 


10 


13 


12 


17 




43 


11 


13 


7 






20 


12 


13 


32 






45 


13 


13 


12 


17 




43 


14 


13 


7 






20 


15 


13 


32 






45 



Table 5 : 15 Yr Quota Calculations Using the Data in Table 4 

and Assumptions in Example 1 . 



Planning 

Period 

0 



School 

1 



24 I 



30 ] 



OD 

nn 



Tour 

3 





0 


0 


0 


30 1 




30 


30 


30 



0 


6 


30 


30 


24 


0 



0 


0 


6 


0 


30 


24 



0 


0 


0 


24 


0 


30 





0 


0 


0 


° 1 




30 


24 


0 



0 


0 


0 


0 


30 


24 




0 


0 


0 


24 


0 


30 




0 


0 


0 


30 


24 


0 



Total 

Billets 



90 



54 



54 



54 



54 



54 



54 



54 



Excess 

Inventory 

0 



36 



6 



0 



0 



0 



0 



0 



Figure 3 : Simple Example Assuming all Billets Removed from People 

at the End of Their Tours. 
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billets are "removed." In reality these billets usually have their P-code 
removed so they are not counted in the total which is used to plan the future 
educational input. The billet still exists usually, but in an uncoded form. 
This is an important point as we shall see. 

The important question is, of the 36 billets which are removed, how 
many are filled with people currently in their first, second or third year in 
the billet? In order for the flows in Figure 2 to be correct we have assumed 
that of the 36 billets removed, 30 had people in their last year in the billet , 
and the other 6 had people in their second year in the billet . This assumption 
leads to the second line in Figure 3. In tour year one there is no spare 
inventory and 30 people in P-coded billets. In tour year two there are 6 now 
in uncoded billets and 24 in P-coded billets, and in tour year three there are 
30 now in uncoded billets. As these move through the system in successive 
planning periods the same school input is generated as in Figure 2. 

The removal of P-codes from billets is done independently of the length 
of time that the person has been in the billet. If 36 billets have their P- 
code removed, it is more realistic to assume that these are equally spread 
over the tour years. Therefore, of the 36 billets removed we assume 12 are 
currently filled by people in their first year, 12 by people in their second 
year, and 12 by people in their third year. This assumption leads to the 
numbers in Figure 4. 

In period 1 there are 12 in "excess inventory" and 18 in P-coded billets 
in tour years 1 through 3. In period 2, 18 P-coded billets become vacant. 

If these are to be filled with new school input then we need 18 as the quota 
in period 1. Following this argument the reader can easily see that the 
quotas obtained in Figure 4 are all equal to the stationary value of 18. 
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Figure 4 : Simple Example Assuming Billets Removed Uniformly 

Across Tour Years. 
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Thus, by making the more realistic assumption that the billets removed are 
distributed uniformly over the tour years we eliminate the cyclic nature of 
the quota and produce a stationary school input. 

The reader will notice that the total excess inventory in Figure 4 
(62) exceeds that in Figure 3 (42). But the question is, how much of this 
increase is unavoidable? Suppose that the billets which have their P-code 
removed are distributed uniformly among tour years. It might be possible to 
use some of the excess inventory in P-coded billets as vacancies occur in 
future years. Let us look at the 12 excess in tour year one, period 1. Can 
these be used to offset new school input in period 1 by being moved to P-coded 
billets in period 2? These 12 people have just started a tour and are in their 
first year. To use them in a P-coded billet would mean transferring them 
after no more than one year in their current billet. They would also be one 
year off in their career path if they were to be kept in a P-coded billet for 
a full 3-year tour. Such movements can be made, but the costs can be high, 
both in the dollar cost of transfer and in morale and efficiency costs asso- 
ciated with broken tours. 

The real system is, of course, more complex, with a two-year school 
period, a 2-3 year intermediate tour between school and P-coded billet, con- 
tinuation rates sometimes less than 1.0 and multiple P-coded tours in a career. 
However, the fundamental problems shown in Figures 3 and 4 still apply. 

Although there appears to be considerable unused P-coded inventory, much more 
of it is unusable than is assumed in the simple model leading to Figure 2. 
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Example 4 : Smoothing the Quotas in Example 2 . 



We return now to example 2 with the data for the 8510 URL P-code. 

Using the data in Table 4 but assuming billets are uniformly distributed over 
the tour years as in example 2, we obtain the quotas in Table 6. The mathe- 
matical description of this smoothing is described in Appendix (iii) together 
with an APL computer program and an example. 



The reader can see immediately how the quotas have been smoothed and 
reach steady state (12, 18, 6, 0) in six years. 
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13 


18 


6 
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Table 6 : 15 Yr Quota Calculations Using Data in Table 4 and 

Assumption in Example 3 . 
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4. Further Smoothing Arguments . 



Table 6 in section 2 illustrates the quotas for the 8510P URL code 
when the smoothing model is used. Although there is much less fluctuation 
between years than is seen in Table 5, a comparison of the quotas in Table 6 
with past school inputs illustrates dramatically the effect of a large billet 
reduction made in one time period with no thought given to its effect on the 
system. 

Past billet requirements have fluctuated, and growth factors have been 
applied to them to forecast future billet requirements. Recall that there is 
typically a five year period between school entrance and P-coded billet entrance. 
Therefore, the quota for next year (FY76) is aimed at meeting requirements in 
FY 81. Again using the 8510P code as an example, growth rates as high as 
23% per year were applied to current billet levels to project ahead five years. 

At the time these growth rates were considered reasonable, given the results 
of a del phi- technique method used by the Navy to determine them. Today the 
growth rates are considered to be zero, and the billet base to which the growth 
rate applied has dropped from about 280 to below 200 (see Figure 7, page 35). 

Table 7 shows past school inputs (at Monterey) for the 8510P code from 
FY 70 through FY 75, together with the forecasted future quotas for FY 76 
through FY 80 from Table 6. Percentage changes from year to year are given 
also. This data is plotted in Figures 5 and 6. 

The curves in Figures 5 and 6 clearly demonstrate what typically happens 
to future quota predictions when large billet reductions occur in a single 
planning period. We see that the system "overreacts" causing severe percentage 
changes from year to year as seen in Figure 6. This phenomenon has been widely 
observed in industrial production processes which try to adjust to changing 
demand. There is a large literature in "Production Smoothing" where the aim 
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Fiscal 


USN 


Percent 


Year 


Input/Quota 


Change 


1970 


96 





71 


82 


-14.6 



72 


74 


- 9.8 






73 


55 


-25.7 






74 


45 


-18.2 


Paj 


;t 




* 








75 


38 


-15.6 


Inputs 


76 


11 


-71.1 


Projected 


77 


22 


+100.0 


Quota 


78 


20 


- 9.1 






79 


25 


+25.0 






80 


29 


+26.0 







Table 7 : Past Inputs and Future Quotas for the 8510P Code . 
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Figure 5 : Past Input and Future Quotas for the 8510P Code . 
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Figure 6 . Percent Change in Past Inputs and Future 8510P Quotas . 



is to even out production so that serious disruptions in the production process 
are avoided. 

In order to use the mathematical concepts of production smoothing to 
an institute of higher education one must know the "costs" of certain disrup- 
tions in the institution. These costs are in terms of factors such as quality 
of education, morale, stability and long-term goals of the faculty, and the 
ability of the faculty to make short term (1-3 years) commitments to research. 
Unlike machinery, which can be shut down and re-started for reasonably well 
predicted costs, we currently cannot measure quantitatively shut-down and 
start-up costs in a graduate school faculty. Any mathematical formulation 
at this stage using production function smoothing in graduate education would 
be sterile. Rather, we proceed with ad hoc approach based on the quota results 
of the models in sections 1 and 2. 

Looking again at Figure 5 we see that if we ignore for the moment 
period FY 76, the quota continues to decrease, "bottoms-out" at 20 in FY 78, 
then gradually increases to a steady-state value of about 36. If we draw a 
straight line from the point for FY 75 to the point for FY 77 we obtain a quota 
of 30 for FY 76 which is about 83% of the steady state value. The effect of 
a quota of 30 in FY 76 on the yearly percentage change curve in Figure 6 is 
shown with the dashed line. Yearly percentage decreases remain in the region 
of 15%-25% through FY 77 (2 years), and then decrease, turn positive for a 
few years, and eventually go to zero. 

What we have done is to use the smoothed transient model in section 2 
to determine quotas for a future 5 year period (or longer if necessary), plot 
the quota together with past school inputs, and use "eye-ball" smoothing, to 
eliminate wild oscillations in the year percentage change curve. In our 
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example a change to the quota model figures is necessary only in the first 
year (FY76). This will usually be the case due to the dampening effects of 
numerous stochastic elements in the real system. Such smoothing is not 
possible with a model which determines the quota for only one future period. 

It is essential for good planning that not only the immediate year's quota 
is determined, but also the quotas for at least a five year period so that 
magnitudes and direction can influence the current quota. It is almost cer- 
tain that the current billet requirements will not stay constant over the next 
5 years. But this does not mean that we should not try to calculate or be 
influenced by future quotas based on these current requirements. When deci- 
sions on factors such as academic tenure and departmental research programs 
affect academic planning over very long periods it is vital that we consider 
future projected quotas. 



26 



5. Inventory/Bi net Ratios . 



A common indicator of the "health" of a given P-code community is the 
ratio of total inventory to billets for a given rank. The purpose of this 
section is to point out some of the problems with this indicator and how cer- 
tain simple modifications might improve its usefulness. 

Before considering the real system let us look back to the simple 
transient problem in Figure 2. Let K(t) be the ratio of inventory to billets 
in period t. In our simple example K should be 1.0, which is the case 
for t = 0. Using the simple ratio of total inyentory to billets we obtain 

K(0) = 1.0, K(l) = 1.7, K(2) = 1.1, K(t) = 1.0, t ^ 3. 

Consider now the arguments leading to Figure 4. We should really be interested 
in the ratio of usable inventory to billets. Let us assume that any inventory 
in Tour year 1 is usable (even though there are not enough billets available 
in which to use them). Then the ratio would be 

K(0) = 1.0, K(l)=^=1.2, K(l) = 1.0, t ^ 2. 

These ratios are much more indicative of the availabilities of P-code inventory 
The first point then is that usable inventory should replace total 
inventory in the numerator of the ratio. But why are we interested in such a 
ratio? Usually it is used as an indicator of future graduate education require 
ments. In this case we should not be interested in the current ratio since any 
new input to graduate school cannot be used in billets for 5 years. A more 
meaningful ratio is 

Usable inventory predicted in 5 years 
Billets predicted in 5 years 
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for the given P-code, designator and rank. In a perfect system this ratio 
should be 1.0, and any deviation from it would indicate that either the real 
system cannot attain this ideal figure because of possibly unavoidable factors, 
or that temporary deficits or excesses exist. 

In order to use this ratio we must be able to determine the usable 
inventory predicted in a given future year. But the effort needed to do this 
is the same effort needed to calculate the quota; in fact the ratio can easily 
be printed out using either QUOTA or SMQUOTA (see Appendix, sections (ii) and 
(iii)). The main advantage of using such a ratio is, of course, that it be 
simple to calculate. 

The simplest modification to the ratio currently used is simply to look 
at the total inventory predicted in 5 years to billets predicted in 5 years. 
Although much simpler to calculate than the ratio above its usefulness in a 
transient situation following a large reduction in billets is not clear. 

Table 8 shows the ratios of predicted usable inventories to billets 
and predicted total inventory SMQUOTA (see Appendix (ii)). Usable inventory 
is assumed to be those in a billet or in their 14th year of service and not 
in a billet. It is assumed that a person is a LCDR only when he has years of 
service between 11 and 16 inclusive. 

The only ratios in this table which can be affected by quota input for 
LCDR are for periods exceeding 5. The first column shows there is no excess 
usable inventory using the quota in Table 6. The second column shows how 
total inventory to billet ratio decreases from 2.40 to 1.97. The figures in 
parenthesis show how the ratios would change if the FY 76 quota were 15 LT's 
and 15 LCDR's, given a total of 30 as shown by the dotted lines in Figures 5 
and 6 of section 4. 
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Planning 

Period 


Usable Inventory/ 
Billets 


Total Inventory/ 
Billets 


1 


1.30 




3.69 




2 


1.33 




3.67 




3 


1.43 




3.54 




4 


1.25 




3.21 




5 


1.15 




2.87 




6 


1.0 


(1.11) 


2.40 


(2.52) 


7 


1.0 


(1.0) 


2.13 


(2.25) 


8 


1.0 


(1.0) 


1.97 


(2.09) 


9 


1.0 


(1.0) 


1.97 


(1.97) 


10 


1.0 


(1.0) 


1.97 


(1.97) 



Table 8 : Predicted LCDR Inventory/Bi llet Ratios Using the Data in Table 4 . 



A problem with the use of ratios is knowing what a reasonable ratio 
should be. The steady state model in section 1 can be used to determine the 
long range ratio for a given policy. Recall that is the number of P-coded 
billets in rank i (for the given P-code and designator), T^. is the tour 
length, and 3^- is the fraction who get to serve a tour. Let L^. be the 
"lifetime," or total expected time an officer spends in rank i. Then in steady 
sta te 




i 

I 

j=l 



X. 

J 



1 

n 

k=j 



Oti. 



L/Bi, 



i = 1,2, 3, 4. 



(5) 



gives the ratio of total inventory to billets. The quotas {X^} are given 
in equations (l)-(4) of section 1. Using equation (5) and the data in Table 2 
of section 1 we find that 
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kJ = 1.97. 



This agrees with the long-run ratio (after at least 8 years) found in Table 8 
using the smoothed transient model. 

Now let kV be the steady state usable inventory to billet ratio in 
rank i. Then one can show that 



KV = 1 + 






1 + a2(l-32)X2/B2 






1 - ^ 



f®2 '1 

j- - J 



if X2 > 0 



if X2 = 0, 



1 + a2(l-33)X2/B2 



- 



1 - 



if 



VlVl ■ V2M 



and 



1 + a^(l“3^)X^/B^ 



K. = 



i-f 



'®4 

TT • - “2l^2^3’‘2 • “ 31 ^ 3 ’' 3 

4 



For the data in Table 2 we find that 



X3>0 

X3 - 0, 



if X4 > 0 
if X4 = 0. 



Kj =1, K2 = 1, K3 = 1, and kJ = 2.95. 

Note that agrees with the ratio in Table 8 after 6 years. 

It is easy to add these ratios Kj" and as outputs to the steady 
state quota model. 
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6. Determination of Requirements . 



The models discussed in this paper produce quotas to meet given P-coded 
billet requirements. The only exception is the smoothing discussed in section 
4 and illustrated by the dotted lines on Figures 5 and 6. Thus it is crucial 
that accurate estimates of future requirements can be made. The 5 year lag time 
from school entry to billet entry adds to the forecasting problem. It is not 
the purpose of this paper to discuss in detail the determination of graduate 
education requirements, but some observations are necessary in order to see the 
quota model in perspective. 

There is no doubt that precise estimates of P-coded billets five years 
in the future cannot be made. The uncertainties in the system introduce large 
variances which can easily be seen in past attempts at forecasting (see Figure 
7). However, some reasons for the variance can be discovered and to some 
extent these can be controlled. For example, all past attempts at forecasting 
have assumed a constant change in billets from year to year (i.e. straight- 
line projection). Extrapolating a growth rate for the immediate year out to 
five years often leads to unreasonable, if not unbelievable, numbers of billets. 

No other models for forecasting requirements have been used. Even now, because 
past forecasts have been so much in error, the method of forecasting future 
billets is to assume no change from current billets over the next five years. 

One possible improvement over this might be to take the current billets in a 
given specialty, multiplied by the forecasted future officer strength (from the 
five year defense plan for example), and divided by the current officer strength. 
This calculation would help to correct for fluctuations in total officer strength. 
These numbers could be further refined by using factors which indicate how a 
given specialty is changing in the Navy over the years. 
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An important point to remember in determining requirements is that it 
should be done independently of the current ability to fill the requirement . 

It is easy to argue for an increase in billets if an excess inventory of 
educated people are available. It is even easier (and currently more convinc- 
ing) to argue that if not enough inventory exists (or will exist) in a field, 
then the billet requirements must be cut because quotas cannot be filled. Such 
arguments lead to neat bookkeeping but avoid the difficult real problems. 

One could argue that no fixed actual requirement exists, and so much 
effort is spent trying to determine non-existing numbers. Recently many billets 
were removed because it was felt that graduate education was not essential to 
the filling of the billet. Such arguments could probably be made on even more 
billets, and the results are probably more a result of the relative obstinacies 
of the two sides, billet removers, or billet keepers, than on any real require- 
ment. At a lower level, if a man has to be able to read numbers off a chart 
and write them down, then clearly it is essential that he can read or write. 

But as educational level increases beyond repetitive trade-type skills, the 
minimum level of education requirements becomes quite fuzzy. This is especially 
so when one realizes that in graduate education the emphasis should be primarily 
on ways of thinking about complex systems , and not simply on learning skills 
at a more advanced level . Not to appreciate this vital difference is to miss 
the point of graduate education. 

The point of this argument is to show that it might be possible to con- 
tinue a given billet without a P-code with little expectancy for improvements 
or changes. Alternatively, it might be P-coded because it is planned that the 
billet in future years should change from its present scope to one requiring 
graduate education. Thus both viewpoints would be correct for differing 
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objectives concerning the billet. It is the reason for the billet's existence 
which should be analyzed in deciding whether it should be P-coded. If agree- 
ment cannot be reached on the reason's for the billet, then the requirement 
for P-coding remains uncertain. 

In the past it has been common to ask the person currently in a billet 
as to whether or not it should be P-coded. But how can this person be expected 
to keep separate the current work he has to do in the billet, with what the 
objectives for the billet are five years in the future? He may currently feel 
that most of his time is spent in "fire-fighting" mode for which he does not 
need an advanced degree. Two points should be made here. Firstly, if he 
does not have an advanced degree, he cannot possibly know if it would help him 
or not; and secondly, what is planned for the billet in five years may be very 
different from what is currently being done. In short, by asking current 
billet holders to ascertain requirements tends to continue past policy mistakes 
into the future. Perhaps a more healthy approach is to have a separate body look 
at the objectives of each community and their billets, and determine from these 
whether or not a graduate education is desirable. 

Finally, increasing the complexity of the quota model to better imitate 
real-life personnel movements is wasted effort when the results are so sensi- 
tive to the unknown billet requirements. When uncertainty exists it is much 
better to aggregate where possible to take advantage of the "law of large 
numbers." Variances in forecasts with aggregation tend to be smaller than in 
forecasts without aggregation. The tendency toward even more finely divided 
specialty codes will only make quota determination even more difficult, with 
even less chance of having a good match of people to billets in five years. It 
is this last match that is important, but only if the specialty code on the 
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billet accurately describes the requirement for the billet. A move toward 
more general coding, with substitution among P-codes, would lead to a more 
flexible system, would indicate a greater understanding of graduate education, 
and would probably have beneficial psychological effects in job satisfaction. 
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Percent of Officer End Strength. 




Figure 7 : Actual and Forecasted 8510P and D Code Billets . 
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APPENDIX 



Terminology of the APL programming language is used frequently in this 
appendix. Readers unfamiliar with this terminology are referred to Katzan 
[1], or Gilman and Rose [2]. 

( i ) The Steady State Model . 

The steady state flowrates (officers per year) given by equations (l)-(4) 
in section 1 are calculated by an APL function called SSQUOTA, which is listed 
in Figure Al. 

Syntax : SSQUOTA is monadic function which takes a vector as its right-hand 

argument. The elements of the vector are the numerical parts of the P-codes 
for which one would like to calculate the steady state quota. For example, 

SSQUOTA 9210 8530 9410 8110 

would result in calculations for computer management, computer science, financial 
management and aeronautical engineering. 

Global Variables : The single subscripted notation in section 1 must now be 

double subscripted. So we now have , T|^^. , and Y|^^- . where 

k indexes a particular P-code/designator. For simplicity in what follows it 
is assumed that only the P-code is considered. However, with a slightly 
different interpretation of k the model can be used for any P-code/designator 
combination. 

The global variables required by SSQUOTA are: 

th 

PCV - a vector of the numerical parts of all relevant P-codes, the k — element 
being the P-code. Thus PCV has an many elements as there are P-codes. 
Let this be m. 

■h h 

B - an m X 4 matrix of billets with (k,i) — element . 



Al 
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- an 
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matrix 


of 


tour lengths with (k,i)— 


element T|^^. . 




CF 


- an 
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matrix 


of 


continuation fractions with 


(k-i)— element 


“ki* 


UF 


- an 
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4 


matrix 


of 


utilization fractions with 


(k-i)— element 


^ki- 


RF 


- an 
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matrix 


of 


reutilization fractions with 


(k-i)— element 


^ki 



Function Description : Line 1 checks that all codes in the right-hand argument 

of SSQUOTA are valid (i.e. they are contained in PcV). Lines 2-6 calculate the 
four quota numbers ^ corresponding to the right-hand 

argument of the function. Note that these calculations are made simultaneously 
for all k. There is no looping through lines 2-6. Lines 7-10 format and 
print the output for each k. Thus the program loops using line 11. The 
formatted output uses the APL+ formatting function hFMT. 

Line 12 uses a function AYN (answer yes or no) to ask if the used would 
like to make changes in the data and recalculate. A "w" answer terminates 
the function. A "YES" answer results in a question asking which P-codes would 
the user like to investigate further (line 13). Lines 15 through 22 allow the 
user to input new data. When this is complete the function returns to line 3 
and repeats the calculations. 

A number of error checking devices have been inserted to prompt the 
user of errors in input. The variable LF which appears results in a "line feed" 
to make the terminal input/output easier to read. 

Output : The output of the calculations in SSQUOTA is a table similar to Tables 

1-3 in section 1. An example of the input/output is given in Figure A. 3, with 
the global variable values displayed in Figure A. 2. In this example m = 2. 
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Figure Al . Program Listing of the Steady State Quota Model . 
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Figure A2 . Global Variable Values Used in Example. 
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Figure A3 . Sample Input/Output Using Steady State Quota Model . 



( i i ) The Basic Transient Model. 



Consider the P-code and designator as fixed, and let: 

Ii,(t) = inventory in period t with k years of service, k = 1,2,..., 30, 

^ t = 0,1,2,... 

Period t = 0 is the current period. 

c. = fraction of those with k years of service who remain to have (k+1) 
years of service, k = 1,2,..., 30. 

S. (t) = additions to inventory in period t with k years of service from 
students in an earlier period. 

d = "delay" from entering school to entering the P-coded billet, (d is 
assumed to be at least 4 years.) 

= year of service when rank i P-coded billets entered, i = 1,2, 3, 4. 

u^. = last year of service for rank i P-coded billets, i = 1,2, 3, 4. 

= number of P-coded billets in rank i. 

q^(t) = school quota in period t for rank i, i = 1,2, 3, 4. 

w^. = tour length of rank i billets 

= u . - 5.. + 1 . 

1 1 

We assume that anyone entering a P-coded billet stays in the billet for the 
full tour length. 

The inventories in the planning period t = 0, {I|^(0)} are given. 

Then future inventories are given by: 

I,^^^(t+1) = C|^I|,(t) + S|^^^(t+1), k = 0,1,2,. ..,29, (Al) 

t S: 0, 

where lQ(t) = 0. We define the "legacy" of past inventories to be 
Thus 

I^^j(t+1) = c^I,^(t). (A2) 
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The {S|^(t)} are given for t = 1,2,3, since these are from students 
currently enrolled in graduate school. The {S|^(t)}, t S: 4 will determine 
future school input quotas. Consider future period t + d for fixed t > 0. 
From our assumptions. 



u . 

1 

Total Inventory Legacy in Rank i = '[ I'(t+d). 

k=Jl. 

The addition to inventory in period (t+d) with k years of service is 

"i 

S|^(r+d) = Max [0, B. - I’(t+d)], k = (A3) 



= 0 



otherwise. 



Thus the quota for rank i in period t is 

l.-l 

1 

q.(t) = S. (t+d)/ n c 

^ 1 k=Jl.-d 

1 



i = 1,2, 3, 4 
t ^ 1. 



(A4) 



The procedure starts with {I|^(0)} and {S|^(t)} t = 1,2,3. Equation 
(Al) is used to calculate {I|^(t)}, t = l,2,...,d. 

Equation (A2) is then used to calculate I||,(d+1) and these are used in 
(A3) to find S|^(l+d). These are converted to quotas q^(l), i = 1,2, 3, 4, 
using (A4). The values of S|^(l+d) are now used in (Al) to give I|^(d+1), 
which in turn give I|^(d+2) using (A2). Use of these in (A3) give Sj^(2+d), 
which are used in (A4) to give q^(2). This procedure continues until all 
quotas are determined for the planning period. 

The basic transient model uses an APL function called QUOTA which is 
listed in Figure A4. 
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Figure A4 . Program Listing of the Basic Transient Model. 



Syntax : QUOTA is a dyadic function taking scalars for both left and right 

hand arguments. The left argument is the P-code number and the right argument 
is the planning period in years. For example 



9211 QUOTA 12 



will calculate the quotas for P-code 9211 for the next 12 years. 



Global Variables : The global variables required by QUOTA are: 



PCV - P-code vector (see (i) above). 

INVm - an m x 30 matrix of current (1974) inventories, where (1/7774)^^ is 

the number currently in service with P-code i and j years of service. 

B - an m X 4 billet matrix (see (i) above). 

D - a scalar giving the delay between school entrance and billet entrance 

(lead time). 



BW - an m X 4 X 2 three dimensional array. Element gives the 

lowest years of Service for a billet with P-code i and rank j (i..e 
i- for the particular P-code). (BW)..„ gives the highest years of 

«J I O ^ 

service for a billet with P-code i and rank j (i.e. u- for the 

J 

particular P-code). 



OR 



STUD 



an m X 30 matrix of continuation rates. [OR)., is the fraction of 

’ O 

officers with P-code i with j years of service who stay in to have 
(j+1) years of service. 



an m X 30 X 4 three dimensional array, where 
additions to inventory in year k in P-code i 
service from currently enrolled students (for 
zero) . 



(5ri/Z?).j|^ is the 
with j years of 
k = 4 all elements are 



Function Description : Lines 1-3 set up various arrays to be used in the func- 
tion. Line 2-8 essentially calculate the S|^(t+d) in (A3). Line 9 calculates 
the new inventory using (Al) and line 11 calculate the quota using (A4). Lines 
12 and 13 format and print the output. 
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Output : The output is a table with 6 columns. Column 1 gives the planning 

year, columns 2-5 give the four quota numbers and column 6 the total quota. 
There is a row for each planning period. An example is shown in Figure A6, 
and the values of the global variables are given in Figure A5. 

( i i i ) The Smoothed Transient Model . 

Table 5 in section 2 shows a sample output of the basic transient 
model. Clearly this model leads to undesirable cycles in the quota. To 
smooth out these cycles we modify the basic model. 

The same notation as section (ii) is used. In steady state, the number 
of billets per year of rank i (for the given P-code and designator) which 
become vacant is B^./w^. 

We now modify (A3) to 



The remaining equations stay the same. The underlying assumption which leads 
to (A3a) is discussed in section 2 and is not repeated here. 

The smoothed transient model uses an APL function called SMQUOTA which 
is listed in Figure A7. 

Syntax : SMQUOTA has the same syntax as QUOTA. (See (ii).) 

Global Variables : SMQUOTA uses the same global variables as QUOTA. 

Function Description : Essentially the only difference between SMQUOTA and QUOTA 

is in line 8, which now uses equation (A3a) in place of (A3). 

Output : A sample output is shown in Figure A8 using the values of the global 
variables in Figure A5. 



(A3a) S|^(t+d) = Max[0,(B^. 
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Figure A5 . Global Variable Values Used in Example . 
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Sample Input/Output Using Basic Transient Model 
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